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Abstract 

We formulate jV=l super Yang-Mills theory in 3+1 dimensions on a two 
dimensional transverse lattice using supersymmetric discrete light cone quanti- 
zation in the large- N c limit. This formulation is free of fermion species doubling. 
We are able to preserve one supersymmetry. We find a rich, non-trivial behavior 
of the mass spectrum as a function of the coupling g^/JSic, and see some sort of 
"transition" in the structure of a bound state as we go from the weak coupling 
to the strong coupling. Using a toy model we give an interpretation of the rich 
behavior of the mass spectrum. We present the mass spectrum as a function of 
the winding number for those states whose color flux winds all the way around 
in one of the transverse directions. We use two fits to the mass spectrum and 
the one that has a string theory justification appears preferable. For those states 
whose color flux is localized we present an extrapolated value for m 2 for some low 
energy bound states in the limit where the numerical resolution goes to infinity. 



1 Introduction 



In the past years, there has been a tremendous amount of progress in the analytical 
understanding of supersymmetric theories due to the discovery of AdS/CFT correspon- 
dence [1], using "orbifolding" [2], "orientifolding" [3, 4], or others. Recently Armoni, 
Shifman and Veneziano have shown that in the large N c limit a non-supersymmetric 
gauge theory with a Dirac fermion in the antisymmetric tensor representation is equiv- 
alent, both perturbatively and nonperturbatively, to Af=l super Yang-Mills (SYM) 
theory in its bosonic sector [3]. Since for N c = 3 the non-supersymmetric gauge theory 
is just one- flavor QCD, even though we have to keep in mind 1/N C corrections, knowing 
the nonperturbative properties of Af—1 SYM is of great importance from phenomeno- 
logical viewpoint as well. Accordingly, one cannot stress enough the importance of 
being able to perform some nonperturbative numerical calculations for Af=l SYM to 
test the predictions made by Armoni, Shifman and Veneziano. However, this is not an 
easy task by any means. This is because the supersymmetry (SUSY) transformation, 
which is an extension of the Poincare transformation, is broken on a lattice due to the 
lack of continuous spatial translational symmetry and because Leibniz rule does not 
hold on a lattice. Therefore, new attempts to put SYM on the lattice are interesting 
and eagerly awaited. 

There have been some promising approaches along this direction [5, 6, 7] that par- 
tially preserve SUSY. However, there still appear to be some remaining issues before 
the method introduced by Cohen, Kaplan, Katz and Unsal becomes a practical com- 
putational approach [5, 8]. The approach by Sugino [6, 7] has encountered unwanted 
surplus modes in four dimensions (in Euclidean space) [7] and is still waiting numerical 
simulations 1 . For other recent progress in an effort to realize SUSY on a lattice, see for 
example Ref. [11, 12]. 

Recently the authors proposed another method to put SYM on a lattice [13]. 
This approach is in fact not a new idea, rather it is motivated by the idea of the 
"(de) construction" [14] and is a mixture of the two existing ideas; the transverse lat- 
tice [15, 16] and supersymmetric discrete light cone quantization (SDLCQ) [17]. Our 
first attempt was made for 2+1 dimensional M = 1 SYM with one transverse lattice. 
Here we present a formulation for 3+1 dimensional Af—1 SYM with a two dimensional 
transverse lattice in the large N c limit. 

At each site of the two dimensional lattice, we have one gauge boson and one four- 
component Majorana spinor. Adjacent sites are connected by the link variables. All 
these fields depend only on the light-cone time and spatial coordinates and are 
associated with two site indices, say In the large N c limit, however, it turns 

out that we are allowed to drop the site indices for our calculation. This is in some 
sense the manifestation of the Eguchi-Kawai reduction [18]. However, it is well known 
that the naive Eguchi-Kawai reduction encounters a problem due to the violation of 
one of the assumptions made by Eguchi and Kawai [19]. That assumption is the 
U{l) d symmetry. Since we do not have to assume the U(l) d symmetry to justify our 
reduction of the transverse lattice degrees of freedom, we believe that we do not have to 

1 We should note that the topological field approach to constructing a supersymmetric theory on a 
lattice utilized by Sugino was first discussed by Catterall in Ref. [9] , which investigates theories without 
a gauge symmetry. Very recently, Catterall has proposed a geometrical approach to M=2 SYM on the 
two dimensional lattice in Ref. [10]. 
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introduce quenching [19] or twisted [20] lattices, which were invented to overcome the 
problem associated with the naive Eguchi-Kawai reduction at weak couplings. For more 
complete and detailed discussion for this claim, see appendix A. With this reduction 
of the transverse degrees of freedom, we can regard all the fields as 1+1 dimensional 
objects. That is to say that we have some complicated 1+1 dimensional field theory 
with some highly non-trivial interactions of the fields. Furthermore, since we can always 
work in the frame where we have zero transverse momenta P X ,P 2 = 0, A/"=l SUSY 
algebra in 3+1 dimensions becomes identical to M=2 SUSY algebra in 1+1 dimensions, 
which is sometimes referred to as A/"=(2,2) SUSY in literature, (2,2) for two Q +, s and 
two Q~'s. We are able to maintain one of this underlying N—(2,2) SUSY algebra in 
our formulation, meaning that we are able to preserve one exact SUSY. 

We discretize light-cone momentum p + by imposing the periodic condition on the 
light-cone spatial coordinate x~ . Thus, we have two spatial lattices and one momen- 
tum lattice in our model. Since we are dealing with spatial lattices, one has to be 
concerned about the notorious fermion doubling problem. In fact it is well known that 
the transverse lattice suffers from the doubling problem [21]. However, the authors 
have found that SDLCQ formulation of a transverse lattice model is automatically free 
of the doubling problem [22]. 

There are some aspects of this calculation that are similar to the 2+1 dimensional 
model [13] and there are others that are not. What is not the same is that the super- 
charge Q~ has terms which have different powers of the coupling g' = g\/W c , where 
a — 1,2. To be more precise, Q~ consists of terms proportional to g' and terms pro- 
portional to g' 3 . The different powers of g' give rise to a rich spectrum as one varies g', 
and the wavefunctions depend on g' . This means that it is possible to see wavefunc- 
tions which are almost vanishing at small couplings, but become very large at strong 
couplings, and vice versa. 

One more thing which is different from the previous case is that our Q~ has terms of 
third and fifth order in dynamical fields, while all of the terms in Q~ are of third order 
for 2+1 dimensional case. This leads to a hamiltonian of eighth order in fields, which 
is of higher order than the hamiltonian of sixth order that we get from the standard 
formulation of 3+1 dimensional Af—1 SYM on the two transverse lattice. We admit 
that this is a disadvantage of our formulation in 3+1 dimensions compared to that in 
2+1 dimensions. Nevertheless, we still think that our approach is more advantageous 
since in the SDLCQ formulation we use Q~, not the hamiltonian, and this Q~ is still 
of lower order in fields than the hamiltonian obtained from the standard formulation, 
and since the standard formulation suffers from the fermion doubling problem. 

Similar to the 2+1 dimensional case we are not able to preserve the full supersym- 
metry algebra. We are able to maintain one exact SUSY. This is attributed to the fact 
that when quantizing the dynamical fields we have to make the link variable, which is a 
unitary matrix, a linear complex matrix. One way to compensate for the effects of this 
"linearization" is to make use of the "color-dielectric" formulation of the lattice gauge 
theory [16, 23, 24]. In this formulation we consider smeared degrees of freedom At, 
which are obtained from the original link variable M by averaging M over some finite 
volume, say YJ ai , M. In order for this smeared theory to be equivalent to the original 
one, we must have an effective potential for the M. defined by integrating out M [23] 

eM-Veff(M)} = f VM5(M -J2M)exp[-S canontcal (M)]. 

av 
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However, this V e ff(A4) can be very complicated and performing the path integral above 
is extremely difficult, if not impossible. Thus, one makes some approximations with 
ansatz to determine V e ff. For more detail, we'd refer the reader to the Refs. [16, 23, 24]. 

To constrain the linearized fields, we require the model to exactly conserve one SUSY 
as we did for our 2+1 dimensional calculation. That is, we present a physical Q~ that 
preserves one SUSY. By "physical" we mean a Q~ which transforms one physical state 
into another physical state. We are not able to fully recover SUSY due to the absence 
of a physical Q+- This defect results in a different number of massless states in the 
bosonic and fermionic sectors. However, we do see the mass degeneracy among the 
massive bosonic and fermionic states. The linearization doubles the bosonic degrees of 
freedom, leading to the SUSY breakdown. The partial recovery of SUSY implies that 
we have cured some but not all of the problems associated with the linearization. 

We are numerically able to identify what we call the cyclic states and non-cyclic 
states by examining the properties of the states. The cyclic states are those whose 
color flux winds all the way around in one or two of the transverse directions. For the 
non-cyclic states the color flux is localized in color space. The cyclic bound states have 
a non-trivial spectrum as a function of the winding number. We find that m 2 for the 
cyclic bound states can be fit by either b + c/Wi + d/Wf or b+cWf + d/W 2 , where b, c, d 
are some constants and Wj is the winding number in the x/-direction with I = 1,2. 
It could be interesting to know how the form of the m 2 changes from weak coupling 
to strong coupling however the complicated spectrum for strong couplings puts this 
beyond our reach at the present time. 

The structure of this paper is the following. In Sec. 2 we present a standard for- 
mulation of J\f=l SYM with a two dimensional transverse lattice and derive constraint 
equations on the physical states. We discuss the implications of those in some detail. 
We give SDLCQ formulation of J\f=l SYM in Sec. 3 and show that this formulation 
is free from the doubling problem. The coupling dependence of the mass spectrum is 
discussed in Sec. 4 followed by numerical results for cyclic bound states in Sec. 5 and 
for non-cyclic bound states in Sec. 6. The summary and possible further directions of 
investigation are given in Sec. 7. Appendix A is to show how we justify the reduction 
of transverse degrees of freedom in the large N c limit. Derived in the appendix B is the 
J\f=l SUSY algebra in Majorana representation in the D + 1 dimensional light-cone 
coordinates with D — 1, 2, 3. 

2 Transverse lattice model in 3+1 dimensions 

In this section we present the standard formulation of a transverse lattice model in 
3+1 dimensions for an J\f—1 supersymmetric SU(N C ) theory with adjoint bosons and 
adjoint fermions in the large- N c limit. We work in light-cone coordinates so that 
x ± = (x° ± x 3 )/\^2. The metric is specified by x ± = x T and x 1 = —xj, where I = 1,2. 
Suppose that there are N sites sites in both the transverse directions x l and x 2 with 
lattice spacing a. With each site, say n = we associate one gauge boson field 

A v ^ n {x^) and one four-component Majorana spinor \I/ n (:r M ), where z/, fi = ±. A u ^s and 
\l/ n 's are in the adjoint representation. The adjacent sites, say n and n + ii, where 
%i is a vector of length a in the direction x 1 , are connected by what we call the link 
variables M^(x M ) and M^(a; M ). M^(x M ) stands for a link which goes from the site n 
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to the site (n + ij) and M^(x^) for a link from the site (n + ij) to n. We impose the 
periodic condition on the transverse sites so that A^ sit( . s i I+n = A n , ^/N gites i I+n = \l/ n , 
^N aites i I +n = and Mj sitesii+n = . Under the transverse gauge transformation 
[16] the fields transform as 

^^M:[/t-i[//[/t Mi-^UnMiui^ * n ^f/At/* (i) 

where g is the coupling constant and U n = U n (x^) is & N c x N c unitary matrix. In all 
earlier work on the transverse lattice [16] \l/ n was in the fundamental representation. 
The link variable can be written as 

M^) = exp (iagA n+il/2 j{x»)) , (2) 

where A n j is the transverse component of the gauge potential at site n and as a — > 
we can formally expand Eq. (2) in powers of a as follows: 

M„V) = 1 + iagA n ^) + ^ [i^^M - <? 2 (A,/(^)) 2 ] + 0(a 3 ). (3) 

In the limit a — > 0, with the substitution of the expansion Eq. (3) for M^, we expect 
everything to coincide with its counterpart in continuum (3+l)-dimensional theory. 
The discrete Lagrangian is then given by 

£ = trj - IfJTF^ + ^(D.M^Miy 



^J2( M X +il M^ +ij M^ - 1) + * n ir»D^ n 



+ ^ n Y\M I n ^ n+lI M I J - M^n-iM-i,)}, 

where the trace has been taken with respect to the color indices, F n ^ v = d^A^ — 
], /i, v = ±. We choose Majorana representation where Majorana 
spinors have real component fields and T's are given by 

r+ _ r° + r 3 / o o\ r°-r 3 / ^ \ 

" v/2 V ^2 J ' - v/2 v J ■ 
The covariant derivative _D M is defined by 

= d^ n + ig[A n ^,y n ], 
D^Ml = d^Ml + igAnjMi - igM^A n+il ^ ^ iagF^ + 0(a 2 ), 
(^M^)t = d»M? - igM?A» + igA» n+ii M? ^4 iagF ^ + 0(a 2 ). 

In the limit a — > we recover the standard Lagrangian as expected. Of course the 
form of this Lagrangian is slightly different from that in Ref. [16] since the fermions 
are in the adjoint representation. This Lagrangian is hermitian and invariant under 
the transformation in Eq. (1) as one would expect. 
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The following Euler-Lagrange equations in the light cone gauge, A n _ = 0, are 
constraint equations. 



d 2 _ A' = gJ+ ^ ig[Ar, d_A T ] + djd^Aj + 2g^ R , (4) 



where 



+ + tf B =_^^, (5) 

Pi = (Ji, /?2 = (73 and V'l^ are the two-component left-moving, right-moving spinors. 

Since these equations only involve the spatial derivative we can solve them for A~ 
and ipLm respectively. Thus the dynamical field degrees of freedom are Mj[, and 

Eq. (4) gives a constraint on physical states \phys), since the zero mode of J + 
acting on any physical state must vanish, 

J+ \phys) = J dx J+(x^)\phys) = for any n = (6) 

This means that the physical states must be color singlet at each site. 

It is straightforward to derive P^ = J dx~T +± , where is the stress-energy 
tensor. We have 



P+ = « 2 E / dx ~^ [4^ d - M n d - M n + ^Rnd-lpRn) , (7) 
n J \ Qj Q J 



p- = a z \ dx'ti 



Aa 2 g 



(8) 



where one should notice that we've kept the non-dynamical fields in the expression for 
P~ to make it look simpler. When one quantizes the dynamical fields, unitarity of 
is lost and becomes an N c x N c complex matrix [16]. One way to compensate for 
the effects of this "linearization" is to make use of the "color dielectric" formulation of 
the lattice gauge theory [16, 23, 24]. We will approach this issue using supersymmetry 
as we've done for the 2+1 dimensional case. 

Having linearized M B , we can expand M l n and ipR n in their Fourier modes as follows; 
at x + = 

1 r°° j- j- 

<„(*") = Y^k dk+ Krs(k + )e- ik+x ~+b^Jk + )e ik+x -), (10) 

where r, s indicate the color indices, ipR n = I U % J, a — 1,2, a} lsr (k + ) creates a link 
variable with momentum k + which carries color r at site n to s at site (n + ii), dl lsr (k + ) 

5 



creates a link with k + which carries color r at site (n + ij) to s at site % and b^ sr creates 
a fermion at the site n which carries color r to s. Quantizing at x + = we have 



[M^ rs (x ),n J Mkl (y )] = 



l -8(x -y ) ——5 rp 5 sq 5 IJ , 

Zi CI a 



"■ rM h 2a?g* ., u a 

6 5 

{Kj,rs( X ~),^kl,p q (y~)} = { U ij,rs( X ~), iu kl,pq(y~)} = ~ S ( X ~ ~ V~) M^<W , (12) 

_ Lb Ch 

where tim, Tty are the conjugate momentum for M, if), respectively, and we wrote out the 
site indices for clarity. Note that we divided 5ik and 5ji by a because 5^ /a — > 5{x l — y 1 ) 
and Sji/a — > 5(x 2 — y 2 ) as a — > 0. Then, one can easily see that these commutation 
relations are satisfied when a's, gPs and fc's satisfy the following: 



a a 



{% S (A: + ), = S ( k+ - P + )77M^. (14) 

with others all being zero. Physical states can be generated by acting on the Fock 
vacuum |0) with these a 7 ^ 's, rf^'s and b aj( 's in such a manner that the constraint Eq. (6) 
is satisfied. 

Before discussing the physical constraint in more detail, let us point out the fact that 
this naive Lagrangian formulation is not free from the fermion species doubling problem, 
while our SDLCQ formulation that we will introduce in the next section actually is [22] . 
Nonetheless, the constraint equation would still be valid since the constraint equation 
(4, 6) was derived from — d + X J C . - = in which we do not have any problematic 
terms responsible for the doubling problem, i.e. the terms which contains the difference 
between fermions at different sites. Therefore, we assume that this physical constraint 
is valid for our SDLCQ formulation in the next section and we will fully utilize it when 
we carry out our numerical calculations. 

With this subtlety in mind, let us complete this section by discussing the physical 
constraint (6) in more detail. The states are all constructed in the large-iV c limit, 
and therefore we need only consider single-trace states. In order for a state to be 
color singlet at each site, each color index has to be contracted at the same site. As an 
example consider a state represented by \phys 1) = d^ rs a^ sr \0), where we've suppressed 
the momentum carried by a 1 ^ and d 1 ^ and we'll do so hereafter unless it's necessary for 
clarity. For this state the color r at site n is carried by a J J to s at site (n + ii) and 
then brought back by d 1 ^ to r at site n. The color r is contracted at site n only and 
the color s at site (n + ij) only. Therefore, this is a physical state satisfying Eq. (6). A 
picture to visualize this case is shown in Fig. la. Diagrammatically, one can say that 
at every point in color space one has to have either no lines or two lines, one of which 
goes into and the other of which comes out of the point, so that the color indices are 
contracted at the same site. 

One also needs to be careful with operator ordering. One can show that the state 
dn,rs a n,stbn,tr\fy * s physical, while the state ^| rs a^| st (i^| tr |0) is unphysical. This statement 
is almost obvious when one recalls what each creation operator does. 
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(a) (b) 

Figure 1: (a)The color charge for the state \phys 1) = d I J rs a I J sr \0). The planes represent 
the color space, a 1 ^ carries color r at site n to s at site n+ij and d J J carries it back to r at 
site n. (b) The color charge for the state \phys 2) = ^(at^-i)^,™ ' ' ' a i t +i J ,te a nt«-l )- 
The lines which intersect a circle represent the color planes at sites. The color goes all 
the way around the transverse lattice. 

We should however note that a true physical state be summed over all the transverse 
sites since we have discrete translational symmetry in the transverse direction. That 
is, for example, the states d[\ rs a[\ sr \0) and d[\ rs ai 2 sr |0) are the same up to a phase 
factor given by exp(iP 2 a). We set the phase factor to one since we take physical states 
to have P 1 = P 2 = 0. The physical state \phys 1) is in fact Y^fli ^ij>s%- jSr |0) with 
the appropriate normalization constant. From a computational point of view this leads 
to a great simplification in the large N c limit. Because as shown in appendix A it turns 
out that in the large N c limit we can drop the site index n from the expression of the 
supercharges and thus can practically set N si t es = 1 for our calculation. This is in some 
sense the manifestation of the Eguchi-Kawai reduction [18]. Eguchi-Kawai reduction 
tells us in the usual lattice theory that the large N c limit allows us to work with only 
one site in each of the space-time directions in Euclidean space. However, the way we 
justify this reduction in our transverse lattice formulation is quite different from the 
way Eguchi and Kawai do in the usual lattice formulation. Therefore, we believe that 
we do not have to introduce quenching [19] or twisted [20] lattices to overcome the 
problem that the naive Eguchi-Kawai reduction comes across at weak couplings [19]. 
We refer the reader to appendix A for more detailed support for this claim. 

Periodic conditions on the fields allow for physical states of the form \phys 2) = 

^n^N^-iii^ru - ■ ■ a i+ii,ts a i\sr\°)^ The color for this state is carried around the 
transverse lattice, as shown in Fig. lb. We will refer to these states as cyclic states. 
The states where the color flux does not go all the way around the transverse lattice we 
will refer to as non-cyclic states. We characterize states by what we call the winding 
number defined by Wi = ni/N sites , where nj = J2n( a n a n ~ d^d^). For N sites = 1, the 
winding number Wj simply gives us the excess number of a 1 ^ over d 1 ^ in a state. We 
use the winding number to classify states since the winding number is a good quantum 
number commuting with Psdlcq as we wn ^ see m ^ ne nex ^ section. In the language 
of the winding number the non-cyclic states are those states with Wi = and cyclic 
states have non-zero Wj. 
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3 SDLCQ of the transverse lattice model 



The transverse lattice formulation of Af = 1 SYM theory in 3+1 dimension presented 
in the previous section has several undesirable features. First and foremost the naive 
Lagrangian suffers from the fermion species doubling problem [22]. Second, the super- 
symmetric structure of the theory is completely hidden. Lastly, the resulting Hamilto- 
nian is 6 th order in the dynamical fields. From the numerical point of view a 6 th order 
interaction makes the theory considerably more difficult to solve. In Ref. [13] we found 
that the (2+l)-dimensional supersymmetric Hamiltonian is only 4 th order making this 
discrete formulation of the theory very different. Unfortunately, it seems this is not the 
case for 3+1 dimensional model. Instead we seem to have supersymmetric Hamiltonian 
of 8 th order in fields. However, since this SDLCQ Hamiltonian is free from the doubling 
problem [22] and since the supercharge Q~, where a — 1, 2, is of 5 th order and it is this 
Q~ that we make use of for our calculations, we think that this SDLCQ formulation 
is still more advantageous than the naive DLCQ formulation. There can, of course, 
be many discrete formulations that correspond to the same continuum theory and it is 
therefore desirable to search for a better one. 

In the spirit of SDLCQ we will attempt a discrete formulation based on the under- 
lying super-algebra of this theory, 

{Q±,Q±} = 2V2P ± 5 a(3 , {Q+,Q«} = 0, (15) 

where a,/3 = 1,2 and the supercharge Q is given by 



Q = Y, dx Jn 



( Qi 


\ 


Qi 




Qi 




K Q~2 


) 



with jn being the supercurrent at the site n = which is a Majorana spinor. 

For the derivation of the super-algebra in Majorana representation Eq. (15), see the 
appendix B. We've set P 1 = with 7 = 1,2 since we're considering the physical states 
only with P 1 = 0. Note that this choice of P 1 has made Eq. (15) coincide with the 
Af=2 super-algebra in 1+1 dimensions also known as A/"=(2,2) super-algebra although 
we are considering J\f=l SYM in 3+1 dimensions. 

In this effort however there are some fundamental limits to how far one can go. As 
we discussed in the previous section the physical states of this theory must conserve the 
color charge at every point on the transverse lattice. Experience with other supersym- 
metric theories indicates that each term in has to be either the product of one M n 
and one tp n or of one Ml and one ip n therefore is unphysical, by which we mean that 
Q+ transforms a physical state into an unphysical one, so that (phys\Q+\phys) = 0. 
While this is not a theorem, it seems very difficult to have any other structure since 
in light cone quantization P + is a kinematic operator and therefore independent of the 
coupling. There appears to be no way to make a physical P + from . We will use P + 
as given in Eq. (7) in what follows. Similarly, we are not able to generally construct 
physical P 1 from and Q~Z. In fact P 1 is unphysical in our formalism, leading to 
(phys\P I \phys) = 0. Formally we will work in the frame where total P 1 is zero, so 
it would appear consistent with this result. We should note, however, that this is not 
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totally satisfying because P 1 = was a choice and a non-zero value is equally valid and 
not consistent with the matrix element. 

Despite these difficulties we find a physical Q~ which gives us Psdlcq Pcont- 
The expression for Q~ is 



i2 3 /V£ / dx-tr 



Mi dl M? + d- M^A - 2g^ Rn 



—i 
2g~a~ 2 

x^M Rn ) a + ^(MlM^M^M^ - l)(fr/W*0c 



I d 3 xl-2gJ + ^-(a 2 ^ R ) a + Fj^tfjo^n) 



where f3i = <7i, /3 2 = 0-3, <7</ + = i<?L4/,cLAj] + d/cLA/ + 2g^ R vjj R , and the last line is 
the continuum form for Q~ in 3+1 dimensions. 

It is tedious but straightforward to check that {Q\,Q\} 7^ {Q^jQ^}"! while both 
{Qi,Qi} and {Q2,Q^,} gi ye the same correct P~ in the limit of a — > 0. In addition, 
one can show that {<5i~> Q2 } 7^ in the discrete form but becomes zero as a — > 0. This 
means that we preserve on/?/ one supersymmetry algebra, say {Qr^Qi - } = P~i i n our 
discrete formalism. We cannot use both Q\ and Q 2 at the same time to construct 
physical states since they do not commute with each other. However, both Q± and 
Q2 separately give us the same mass spectrum when we perform SDLCQ calculations. 
Thus, it is suffice to consider only one of the two and we take Qi for our calculations 
in the following sessions. 

Notice that Q~ above is fifth order and, thus, Psdlcq obtained from it is eighth 
order in fields as we mentioned at the beginning of this section. In fact we find 



SDLCQ 



= « 2 E/ 



{QI,QI} 



2V2 



dxtr< 



9" 



Mi d Ml} + MlU, d Ml_ lr ) + 2« 



x 



d 2 2g 2 a 2 



2 2g 2 a 2 

Ml d M? + Mll H d Ml_ ir ) + 2« 



u'M^-MSDd-'iMlul 



n—ij 



—{ 



n+iz "ii- -in ■ 

)d~ 1 (M 2 Ju 1 n M 1 n M 2 n+ti - u\ +l M 2 n+l Ml +2l Mll ll+ ^ 
+ (ul +l2 M 2 J - M^u^d-\u\MlMl Mx M% 2 - MiUX-nM'-iM-n+iJ 
Hu 2 n+il M 2 n+h - Ml+ i y n+li+t2 )d~\Ml\ t2 M^u l n M 1 n - Ml +li+l2 Mll 2n Ml\ n u l n+ti ) 
H<Ml - My n+n )d-\M 2 n+l MH l2 M^ul - Mll n _ l2 Mll t2 ul_ l2 M 2 n _ l2 )} 



H+il+12 "+12 



■ II 

n+12 



Aa 2 



{{M 2 n+li Ml\ l2 M^u l n 



M 2 ^ M 1 ^ 11 ■ M 2 A 



x0->^M^XU - Mll l2 u l n _ l2 M l n _ l2 M 2 n+n _ l2 ) 
+ {M^u l n M l n M 2 n+ti - ul +i2 M 2 +h Ml +2i2 M 2t +h+i2 ) 
xd-^M^M'JulMi - M 2 n+ii+i2 Ml\ 2i2 Ml\ i2 u l n+i2 ) 
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+ {MH l2 M^u l n M l n - Ml +ii+i2 Ml\ 2h M^u~ n+il/ 



rlt 



U n+i 1 Mn+i! Mn+2i! M n+ h +i 2 , 



xd-^M^ulMlM 1 ^ 
+ {u l n M l n M 2 n+li M l n \ l2 - M^X-iX-iX-ii-to 



-1 

Ik 2 



(M l n Ml +l M l n \ t2 M% - MlMl +l Ml\ n M^f 



One can show that by setting g = and M, AP = 1 this Psdlcq gi yes r i se to a 
dispersion relation 



AT = 



2k+ 



sin 



a/2 



+ 



sin 



a/2 



which is free from the fermion species doubling problem [22]. Furthermore, one can 
check that this Q~ commutes with P + obtained from C; [Q~ , P + ] = 0. Thus, it follows 
that, 

(phys\[Q~,M 2 }\phys) = (phys\[Q~ , 2P + Ps DLCQ \\phys) = (16) 

in our SDLCQ formalism, where M 2 = 2P + P SDLCQ - (P 1 ) 2 - (P 2 ) 2 . The fact that the 
Hamiltonian is the square of a supercharge will guarantee the usual super symmetric 
degeneracy of the massive spectrum, and our numerical solutions will substantiate this. 
Unfortunately one needs a Q + to guarantee the degeneracy of the massless bound states. 

Recalling that we set N s n es = 1 in both transverse directions and that we are in the 
large-iV c limit, we can write Qi as 

Ql = Qu + Ql2 + 



where 



22" 1 /4 a 2 
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71 



dkidk 2 dk 3 5(ki + k 2 — ks) 



x 



-fPd'a 1 + d'W - b^a'd 1 + aW) 
-d 7 W + o 2 W - a £ W + o 2 W) 



k 2 -h ( 

| k 2 + h 
k 1 y / k 2 k 3 

+ h tb_ (a I Wb 2 - a^a 1 + d'W - d^d 1 ) 

+ ( 1 + 1 - p\ {b 2 Wb 2 + o 2 W) 

\ ki k 2 k 3 J 



(17) 



Qr 2 = 



i2- 1 ' i a 2 g roo 



a g f°° 

= — / dk\dk 2 dk 3 5{ki + k 2 — fc 3 ) 

7T JO 



X 



-^( & 2 W + 6 lt 6 lt 6 2 } + _L (6 lt 6 2 6 l + h 2 W) + J_( fe lt 6 2 tfe l + fe l t6 l &2) l ) (lg) 
K3 K2 J 
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i2~ 1/4 a 2 g a 2 g 2 f°° r 
Q 13 = 1= — / dk\dk 2 dk 3 dk 4 dkA 5(k l + k 2 + k 3 + /c 4 - k 5 ) 

y/TT An Jo L 

„ ,\ , (^WaV + d^dWV - b l U 2 d l a 2 a 1 - cWVW) 
\/k\k 2 k 3 k 4 

+ 1 (cPWcPa 1 + tfWWaPd 1 - d 1 ^b l d 2 d 1 a 2 - b^d^d^a^d 2 ) 
\/k 2 k 3 k 4 k^ 

+ „ \ , (rfVfeW + a^^d^d 2 - d 2 ^a l b l d 2 d 1 - a'Wd^d^d 1 ) 
y/k^kik^ki 

+ ; 1 (A^W + trVV^V - a^aW - aV&^ta 2 ) 

+ „ 1 (a^VaV + dWW - a 2 U l a 2 a l b l - dVaW) 
Vhkik 2 h 

+5(ki + k 2 + k 3 — k 4 - k 5 ) 

' ^dWW + a^dVa 2 - d 2 WtaV - a 2 WaW) 



y/kik 2 k 3 k. 
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+ - 1 (fc^tWW + dV&W - b^d^d^dW - d^d^bW) 

Vk 2 k 3 k 4 k 5 

+ ,,},, {a l Wd 2 U 2 a l + dVWd 1 - tWWa 2 - dVaW) 

+ /; ,\ , (a 2 Wt Q y + a 2 WaV - aWVa 1 - aWaV) 
\Jk 4 k^k\k 2 

+ - 1 (WWd* + fo lt rf 2 trf l rf 2 a l _ d l td 2t a lt 6 l d 2 _ & lf d lt d 2 d l a2) l 1 (ig) 

Vk 5 kik 2 k 3 J J 

with A; + = fc, ai = a(fci), a^aa = tr(c4aia2), a^a^a = t^aja^), a^aaaa = ^(05010203(14), 
aVoWa = tr(a{a2 a 3 a 4 a 5)> a^a^aa = t^aJaJjaJc^as), and aWaaa = ti{a\a\aia 2 a 3 ) . 
Qn is the part of which looks exactly like Q~ in 2+1 dimensional model with the 
difference being that here we have two types for each of the bosonic fields a and d. Q\ 2 
is a new piece in 3+1 dimensions and mixes two different types of fermionic fields. Q13 
is also new and composed of fields of fifth order. Note that for small couplings, Qn 
and Q12 dominate over Q 13 , while Q 13 dominates in the strong coupling regime. Notice 
that from this explicit expression for Q± it is clear that the winding number introduced 
in the last section evidently commutes with Q± and, thus, with Psdlcq- Therefore, 
cyclic states do not mix with non-cyclic states. 

It is always important to look for symmetries of Q~ since the symmetries, if any, 
will reduce the amount of the computational efforts considerably. To do this, let us 
consider three cases separately: (i) the intermediate coupling where we have all the 
three pieces together for Q{; (ii) the weak coupling limit where we can ignore <2 i3 ; (iii) 
the strong coupling limit where we consider Q i3 only. For the first case (i) we find two 
Z 2 symmetries, 



• a\j <-»• -a 2 p d\j <-> -d 2 p b\ <-> -b\ p b 2 unchanged, 

• a\j <-> -d^, bfj <-> 

The first symmetry implies that states with the winding numbers, say (W\, W 2 ), are 
equivalent to those with (W 2 , W\) up to the minus sign. On the other hand the second 
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symmetry implies that states with {Wi, W 2 ) are equivalent to those with (—Wi, — W 2 ) 
up to the minus sign. 

In the case of the weak coupling limit (ii), we find two more independent Z 2 sym- 
metries; 

# fl l ^ „2 _ _„2 ^2 ^ _^/2 



The second of these implies, with the help of the second Z 2 symmetry we found in 
the case of (i), the equivalence of states under {Wi, W 2 ) <-> {—Wi, W 2 ) <-> (Wi, —W 2 ). 

In the strong coupling limit (iii), we do not have any other Z 2 symmetries besides 
the two we found in the case of (i). However, it is easy to see that Q 13 commutes with 
b 2 ^b 2 , thus the number of b 2 ^s is a good quantum number as well as the two winding 
numbers. 

It is interesting to see what we can find for each of the three different cases (i), 
(ii) and (iii). However, in this our first attempt to formulate M = (2,2) SYM in 3+1 
dimensions with SDLCQ on a two dimensional transverse lattice, we constrain ourselves 
to consider only the most generic case (i) where we have all the three pieces together 
for Q-. 

Now we are in a position to solve the eigenvalue problem 2P + Pg DLC n\phys) = 
m 2 \phys). We impose the periodicity condition on M^, M^' and u" in the x~ direction 
giving a discrete spectrum for k + , and ignore the zero- mode: 

7T r°° 
k + = - n (n = l,2, 



We impose a cut-off on the total longitudinal momentum P + i.e. P + = irK/L, where 
K is an integer also known as the 'harmonic resolution', which indicates the coarseness 
of our numerical results. For a fixed P + i.e. a fixed K, the number of partons in a 
state is limited up to the maximum, that is K, so that the total number of Fock states 
is finite, and, therefore, we have reduced the infinite dimensional eigenvalue problem 
to a finite dimensional one. 

For this initial study of the transverse lattice we consider resolution up to K = 8 
for non-cyclic {Wi = W 2 = 0) states and up to K = Wj + 6 and K = Wj + 5 for 
states with \Wj\ = 1 and \Wj\ = 2,3,4,5, respectively. We were able to handle these 
calculations with our SDLCQ Mathematica code and C + + code. 



4 Coupling dependence of the mass spectrum 

In this section we will discuss the mass spectrum as a function of g' = g^/W c for 
K = 4,5,6. 

It is instructive to see the dependence of m 2 on the coupling since we have terms 
in Q~ that go like g' and g' 3 . In Fig. 2 we show the entire mass spectrum of non-cyclic 
states in units of g' 2 /na 2 for K = 4, 5, 6 as a function of g' in a log-log plot. In order 
to see the crossings in more detail we show Fig. 2(b), (d), and (f) on a different scale 
from (a), (c) and (e), respectively. We've set 1CT 8 or less to the numerical zero in our 
code. 
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Figure 2: Log-log plots of the mass spectrum m 2 in units of ^2 versus g' = g\fW c with 
K — 4, 5,6 for (a),(c),(e), respectively, (b), (d), and (f) are the same as (a), (c) and 
(e), respectively but on a different scale so that one can see the crossings in more detail. 
10~ 8 or less is the numerical zero in our code. 

As one can see from Fig. 2, there is a rich structure in the mass spectrum as a 
function of g', and the origin of this structure for the case where K = 4 in Fig. 2(a) and 
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(b) is rather easy to understand. We find four types of states; (i) those states which are 
killed by Q13 and whose m 2 in units of ^ are independent of g'; (ii) those states which 

vanish upon the action of Qn + Q 12 and thus whose m 2 in units of ^ go like g' 4 ; (iii) 
those states which survive upon the action of Qn + Q 12 and of <2 13 independently and 
whose m 2 in units of ^ go like (A + Bg' 2 ) 2 , where A, B are some constants; (iv) those 
massless states which become zero upon the action of Q n + Q 12 + Q13. From Fig. 2(a) 
and (b) it is easy to identify one state each for the second and third type because m 2 of 
a state of the second type go like g' A , giving rise to a straight line with a non-zero slope 
for all g' in the log-log plot, while m 2 for the third type is (A + Bg' 2 ) 2 , leading to some 
flat, constant line at small g' and a (inclined) straight line at large g'. We should note 
that for the second kind one should take into account the level crossing. The rest of the 
states clearly fall into either the first kind or the fourth kind. States of the first type 
yield (/-independent m 2 , thus, a flat line in the log-log plot, while states of the fourth 
type are massless represented by the "dots" below the line of log 10 m 2 = —8 since the 
numerical zero is set to 10~ 8 in our code. 

This discussion does not however seem to explain the dependence on g' of the mass 
spectrum with K — 5,6. To get the full understanding of the behavior, we made a toy 
model. In this model we have a 2 x 2 matrix R for the boson sector of Q~ given by 




bi + c l9 ' 2 b 2 + c 2 g' 2 \ 
h + c 3 g' 2 64 + c±g' 2 J ' 



where foj with « = 1,2,3,4 is equal to either or 1 and q is equal to either or 1/47T. 
Here one should notice that we've factored out g' from R or Q~ , and therefore g' 2 from 
P~ . The Q~ for this toy model is thus given by 

<™ = ( £ ) • 

where T stands for the transpose. Thus, the matrix to diagonalize is 

RR T \ 
R T R ) ' 

or equivalently RR T . Among the 2 8 = 256 possible forms for Q~, we found sets of 
parameters that lead to a level crossing, and non-trivial behaviors in the mass spectrum. 
Some of those non-trivial ones look the same as some of those in Fig. 2, while there 
are others which do not look like any of those in Fig. 2. For example see Fig. 3, where 
Fig. 3(a) and (b) are the ones that we can see in the actual spectrum in Fig. 2, while 
3(c) and (d) are not. The sets of parameters we used are given in Table 1. Of course 
there are ones which are seen in Fig. 2, but cannot be found in our toy model. However, 
it is very likely that as we increase the size of the matrix R of our toy model, we would 
be able to identify those not-yet-seen behaviors in our toy model as well. 

Using this toy model, we can study wavefunction dependence on the coupling g'. 
As the simplest example, consider the case of the level crossing shown in Fig. 3(a). In 
this case we can think of a bound state \m 2 ) as a linear combination of two different 
states, 

l™ 2 > = /G/)|l> + %')|2>, 



(Q-/g r Y 
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l°g,„g' log,„g' 



(c) (d) 

Figure 3: Sample spectra obtained from our toy model, (a) and (b) can be seen in the 
actual full spectrum in Fig. 2, while (c) and (d) cannot. 

where f(g') and h(g') are wavef unctions, which depend on g'. |1) is a state of the first 
type of the four we considered above and responsible for the constant behavior of the 
mass spectrum and |2) is a state of the second type responsible for the g /4 -behavior. 
In Fig. 3(a) the higher energy state stays constant for small g', where f(g') ^> h(g'), 
and goes like g' 4 for large g\ where h(g') 3> f(g')- The opposite behavior of the 
wavefunctions is true for the lower energy state. That is, the lower energy state goes like 
g' 4 for small g', where h(g') ^> f(g'), and stays flat for large g', with f(g') ^> h(g'). This 
observation implies that for more general cases a bound state is a linear combination of 
states of the four types associated with ^'-dependent wavefunctions, and it is the non- 
trivial g'-dependence of the wavefunctions that gives rise to such a rich, complicated 
spectrum in Fig. 2. 

We expect that the structure of the mass spectrum as a function of g' will persist 
for the cyclic states and in fact we have numerically confirmed the similar structure for 
them as well. 

Note that since the dominant structure of a bound state changes as one changes g', 
there is some sort of "transition" as one goes from weak coupling to strong coupling. 
It is of great interest to see if the winding number dependence of the mass spectrum 
varies due to this transition. We are not able to identify any states in strong coupling 
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Table 1: Parameter sets used for our toy model to get each of the spectra in Fig. 3. 

regime because of the rich and complicated behavior of the spectrum although we are 
able to find some states in the intermediate region where g' — 1. 

We discuss the mass spectrum of the cyclic states as a function of the winding 
number and the resolution with g' = 1 in more detail in the next section. The discussion 
of the mass spectrum of the non-cyclic states is in the following section. 

5 Numerical results for the cyclic (Wi ^ 0) bound 
states 

In principle we can study the case where both of the winding numbers are non-zero and 
the case where one of them equals zero. However, the size of the Fock basis is much 
larger for the former case than for the latter. This means that we can reach a higher 
resolution for the latter case. Thus, in order to get enough data to analyze for our 
first attempt we restrict ourselves to the case where we set one of the winding numbers 
to zero. Since we have two Z 2 symmetries, (Wi,^) (W2,Wi) and (W±, W 2 ) <-> 
(—Wi, — W 2 ), we can set W 2 = without loss of generality and consider only positive 
Wi when studying the winding number dependence of the bound states. 

As guaranteed by the super-algebra, we find numerically a degeneracy in the mass 
spectrum between massive fermionic and bosonic states. However, this supersymmetry 
is broken for the massless states since we do not preserve the entire set of super sym- 
metry algebra. In this section we only consider the massive bound states, and therefore 
it suffices to consider only bosonic states. 

In Fig. 4(a), (b), (c), and (d) we give plots of m 2 with g' — 1 for four low-energy 
bound states as a function of l/(K — W\) and extrapolate m 2 ^ in the (K — W\) — > 
00 limit using a linear fit b + c/{K — W\) for (a) through (c) and a quadratic fit 
b + c/(K — Wi) + d/(K — Wi) 2 for (d), where b, c, d are fitting parameters. We identify 
a bound state with different K's from the properties of the bound state, such as the 
averaged number of partons of a particular type etc. We present here four bound states 
we could easily identify. The dominate fock component of the bound state in (a) and 
(c) has the form b^a 1 * ■ ■ -a^b 1 ^. For the bound state in (b) the dominant component 
is of the form b^a 1 * ■ ■ -a^b 2 ^ . The bound state in (d) has the dominant component of 

In Fig. 5 we present m 2 ^, obtained in Fig. 4(a), (b), (c), and (d), as a function of 
W\. We show a fit to the data of the form b + cW 2 + d/W 2 in Fig. 5(a) and of the form 
b + c/Wi + d/W 2 in Fig. 5(b). As can be seen, it is difficult to say which fit is better 
from the graphs. The fit of the form b + cW 2 + d/W 2 appears a bit better. 
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Figure 4: Plots of m 2 in units of ^ of low-energy cyclic bound states versus \/{K—W\) 
for 1^1=1 (circle), 2(square), 3(diamond), 4(triangle up), 5(triangle left). Also shown are 
a linear fit for (a), (b), and (c) and a quadratic fit for (d). The coupling g' = g\fW c = 1- 



The use of a fit of the form b + cW 2 + d/W 2 has a string theory justification. In the 
string theory the energy of a string confined in one dimension with a period L is given 
by the sum of its momentum mode and its winding mode, so that E = p2n/L + qTL, 
where p, q are integers and T is the string tension. Now if we consider our cyclic bound 
states as a string confined in the Xi-direction with L = aW±, then it follows that 
m 2 = b + cW 2 + d/W 2 . 

We should however remind the reader that we used a fit of the form b + c/W + d/W 2 
in Ref. [13]. There we argued that the operator has the form Q~ = b + ck± in 2+1 
continuous theory and m 2 ~ (Q~) 2 = b + c/W + d/W 2 with k± ~ 1/L ~ 1/W. This 
behavior is consistent with the unique properties of SYM theories that we have seen 
in previous SDLCQ calculations [25]. We have seen that as we increase K we uncover 
longer bound states that have lower masses. Supersymmetric theories like to have light 
bound states with long strings of gluons. We call these bound states with long strings 
of gluons, stringy bound states. In 3+1 dimensions with two transverse lattices we have 
seen the stringy bound states as well, and we have Q~ = b + ck\ + <iA; 2 , leading to the fit 
of the form b + cjW x + d/W 2 in Fig. 5(b) for ki ~ 1/L ~ 1 /W x and k 2 = 0. Up to the 
numerical resolution we can correctly reach, we can not say for sure which form of m 2 
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Figure 5: Plots of K — > oo limit of m 2 in units of ^2 of low energy cyclic bound states 
versus W 1 with a fit to the data of the form b + cW 2 + (i/W^- 2 in (a) and of the form 
b + c/Wi + d/W 2 in (b). The circles correspond to the bound state in Fig. 4(a), squares 
in 4(b), diamonds in 4(c), and triangles in 4(d). 

describes A/"=(2,2) SYM in 3+1 dimensions. It appears that the form b + cW 2 + d/W 2 
is preferable, suggesting that the cyclic bound states in 3+1 dimensions are more like 
a string with the energy of the form E = p2ir/L + qTL. 

6 Numerical results for the non-cyclic bound states 

(W, = 0) 




1/K 

Figure 6: Plots of m 2 in units of ^ of low-energy non-cyclic bound states against 1/K 
with a linear fit to the data. The coupling g' = g\fN~ c = 1. The circles correspond to 
bound state A, squares to the state B, diamonds to state C 

Let us now discuss numerical results for the non-cyclic bound states. Again we 
follow bound states that we can easily identify from the properties of the bound states. 
In Fig. 6 we show m 2 of three low-energy states in units of ^ as a function of 1/K 
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Table 2: Extrapolated values for m 2 in units of ^ as K — > oo for State A, B, and C 
in Fig. 6 represented by its dominant Fock state. 

with g' = g\fW c = 1. The state A denoted by circles is composed primarily of two 
bosons and two fermions, b^d^a^b 1 ^ . The state B and C denoted by squares and 
diamonds are composed primarily of two fermions, b^b 2 ^ and b^b 1 ^, respectively. We 
show a linear fit to the data and see good conversion as K — > oo for all the three states. 
The extrapolated values for m 2 in the limit of K — > oo are given in Table 2. We also 
find the stringy states for the non-cyclic states. 

Recall that we found in Sec. 4 that a bound state would be a linear combination 
of states of the four types we enumerated in Sec. 4. Hence, it is instructive to see if 
we can identify the three bound states with any of the four types. For K = 4 we can 
identify all the three bound states with those that are killed by Q13 and whose mass in 
units of f^z are independent of g' . However, as K increases, we are not able to classify 
them into any particular type of the four. This is because as we increase K the number 
of states becomes very large and the mass spectrum becomes dense. It is likely that 
these states mix with other nearby states with the same coupling dependence, giving 
rise to small changes in m 2 but still the same general coupling dependence. At this 
time however we are not able to resolve the spectrum in an enough detail to study these 
effects. 

7 Discussion 

We have presented the standard formulation of M=l SYM in 3+1 dimensions with a 
two spatial dimensional transverse lattice. Then we gave the SDLCQ formulation of 
the theory. We found that the standard formulation suffers from a fermion species dou- 
bling problem, while SDLCQ formulation does not. In the frame where the transverse 
momenta equal to zero, M—l SUSY in 3+1 dimensions is equivalent to M=2 SUSY 
in 1+1 dimensions also known as A/"=(2,2) SUSY. We were able to present Q~ which 
has the correct continuum form and yields by the SUSY algebra a discrete form of P~, 
where a — 1,2. This P~ then coincides with its continuum form in the continuum 
limit. Since Qi and Q2 don't commute with each other in our formulation, we are to 
use only one of them to solve the mass eigenvalue problem, preserving one exact SUSY. 

We found that this Q~ consists of terms which are proportional to g' = g\fN~ c an d 
terms which go like g' 3 . This led us to investigate in some detail the g' dependence of 
the mass spectrum. From a simple toy model we concluded that the rich, complicated 
behavior of the mass spectrum with varying g' is due to some non-trivial coupling 
dependence of the wavef unctions. This is also responsible for a "transition" in the 
structure of a bound state when going from weak coupling to strong coupling. Because 
the dominant structure of a bound state changes with changing g' . 

We classified the bound states into two types, the cyclic and non-cyclic as we did in 
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Ref. [13]. The cyclic bound states are those whose color flux goes all the way around in 
one or two of the transverse directions. The bound states whose color flux is localized 
and does not wind around are referred to as the non-cyclic bound states. For each type 
of the bound states, we were able to identify some bound states in the mass spectrum 
for g' — 1 and found the K — > oo limit of m 2 . 

For the cyclic bound states we were able to present m 2 as a function of the winding 
number Wj in the xj direction with I = 1,2. We found two very good fits to the data. 
The first fit b + cW 2 + d/Wj is motivated by the string theory, where the energy has 
the form E = plnjL + qTL, where p, q are some integers, T is the string tension and 
L is the period of the transverse lattice. The other fit b + c/Wj + d/W 2 is motivated 
by the operator structure of Q~. It appeared that b + cW 2 + d/W 2 is preferable, 

For the non-cyclic states as 1/K — > we saw good linear conversion of m 2 of the 
low-energy bound states that we could identify and gave the extrapolated values for 
m 2 . We could identify for K = 4 the bound states with a state whose m 2 in units of 
■^-2 are independent of g' though we were not able to do so for higher K's because of 
the dense, and complicated spectrum. 

In summary, we were able to present a formulation of SYM in 3+1 dimensions with 
one exact SUSY on a two dimensional transverse lattice and find the mass spectrum 
nonperturbatively. There remain however a number of important questions to answer. 
First and foremost it is of great importance to determine the form of m 2 numerically to 
better precision. It is interesting to see what the winding number dependence of m 2 is 
if both of the winding numbers are non-zero. We need to invent a method to resolve the 
dense spectrum at strong couplings. This will help us see if there is any "transition" in 
the form of m 2 as one goes from weak coupling to strong coupling. However, perhaps 
most importantly, as discussed in appendix A we need to know to what extent we've 
resolved the problem caused by the linearization of the link variables that we needed 
to quantize the fields. Knowing this tells us how reliable our numerical results are. 
Because one of our major simplifications in numerical calculation in the large N c limit 
comes about from the reduction of the transverse degrees of freedom whose justification 
relies upon the presence of the quantized fields and the vacuum. Restoration of SUSY 
for massive bound states, which has been broken by the linearlization gives us some 
confidence that our formulation indeed provides some sensible results. However, we 
would still have to clarify the issue to be more certain and confident. To this end, we 
need to compare our numerical results with some well-established theoretical predictions 
and with other numerical results obtained from the usual lattice calculation. Hence, 
it is of importance to apply our formulation to some other supersymmetric theories in 
higher than 1+1 dimensions, for instance, Wess-Zumino model, lattice sigma model, 
and SQED. It appears that the application is relatively straightforward. From more 
practical point of view, a next question to ask is what happens if one includes scalars 
and their superpartners in theory. We did not consider this case in this paper simply 
because this was the first attempt to formulate SYM in 3+1 dimensions with one exact 
SUSY on a two dimensional transverse lattice and, thus, we wanted to consider the 
simplest possible case. However, it is of great interest to consider the question in the 
future. The authors believe that when we are able to answer all those questions, we 
will also be able to test the predictions made by Armoni, Shifman and Veneziano [3, 4]. 
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A Eguchi-Kawai Reduction 

For our numerical calculation we've set N S i tes = 1, in other words, we've dropped the 
site indices. This reduction of the transverse degrees of freedom has brought a great 
amount of simplification in our calculation and needs some detailed justification. Since 
it is only the supercharges that we need to do our calculation, if the supercharges do 
not depend on the site indices in the large N c limit, neither does any quantity that can 
be computed from Q~, for instance m 2 for our case. Therefore, in order to justify the 
reduction of the degrees of freedom for our purposes, it suffices to show the independence 
of Q~ of the site indices in the large iV c limit. In this appendix, in particular, we will 
show that in the large N c limit the leading order terms of the supercharges Q~ with 
keeping all the site indices are the same as those with setting N sites = 1. We should 
note that this sort of arguments about the justification of the reduction on a transverse 
lattice have already been given in literature, for instance see Refs. [16, 23, 24] and our 
arguments below closely parallel those in the Refs. [16, 24]. 

In what follows we only consider Qi, however the same arguments apply equally 
well to Q 2 - For definiteness let us first consider a Fock state denoted by 

5>[. . .di^h^^k^^h^ih) . . .]|o), 

n 

where we've written k + = k, n = is the transverse lattice site, i\ is the vector of 
length a pointing the x 1 direction, a is the lattice spacing, and the dots represent some 
creation operators. When we act on this state with Qi , we get for example from one 
of the terms in Q\, say b 2 ^d}a} = Z) n tr[6^(pi + P2)dl l (pi)al l (p 2 )] on it 

N c J2tr[. . . b 2 }{k 1 + k 2 )b^(h)a l l n (k A ) . . .]|0). 

n 

If we set N sites = 1, then the Fock state now becomes 

tr[. . . d lt (A; 1 )a lt (A; 2 )fo lt (A;3)a lt (A;4) . . .]|0>, 
and upon the action of Q{ we get from b^d^a 1 = ti\b 2 \pi + p 2 )d 1 (p 1 )a 1 (p 2 )} on it 

7V c tr[. . . b 2 \ki + k 2 )b 1 \k 3 )a 1 \k A ) . . .]|0), (20) 

and one more term 

tr[. . . b 2 \ki + k 4 ) . . .]tr[a 1 +(A: 2 )& 1 t(fc3)]|0). (21) 

Notice that the extra term Eq. (21) we get by setting N sites = 1 is down by 1/N C 
compared to the leading order term Eq. (20) and thus we can ignore it in the large N c 
limit. Of course, in the above example, we could and would have gotten many more 
terms depending on what we have in those 'dots' inside the trace of the Fock state we 
considered. However, it is easy to see that our conclusion remains the same; all the 
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extra terms we get by having only one site are down by 1/N C or more powers of 1/N C . 
This all comes down to the fact that we can have only single-traced states in the large 
N c limit. Therefore, we find that the leading order terms of Q~ are the same whether 
we keep track of the site indices or not. Although this proof is for finite K, we suspect 
that the same result would persist at infinite K. 

The way to justify the reduction here should be contrasted to the way exploited by 
Eguchi and Kawai [18]. Eguchi and Kawai showed that in the large N c limit we can 
work with only one lattice site in each of the space-time directions in Euclidean space. 
However, the proof was based on, among others, the assumption that U{l) d symmetry 
is not spontaneously broken, where d is the number of the space-time dimensions. This 
assumption was found to be wrong for d > 2 at weak couplings by the authors of 
Ref.[19]. To resolve this problem, there have been many models proposed, for instance 
quenching [19] and twisted [20] lattice formulations. Here in our formulation, however, 
we believe that we do not have to introduce any of the modified lattice formulation 
since the way we justify the reduction is quite different the way Eguchi and Kawai 
do. Our proof stands on its own feet regardless of us maintaining the U(l) d symmetry 
or not and, therefore, would not suffer from the problem associated with the naive 
Eguchi-Kawai reduction as we go from weak to strong couplings. 

A question, however, remains. That is the question of how well we've managed to 
quantize the fields since all our arguments above rely on the fact that we have the quan- 
tized fields and true vacuum. Put in another way, how good the reduction procedure 
is depends on how good our quantization procedure is. Recall that to quantize, we had 
to "linearize" the unitary link variables, which leads to the breakdown of SUSY. The 
authors of Refs. [16, 23, 24] make use of the "color-dielectric" formulation to resolve 
the problem for non-supersymmetric theories. Although this formulation resolves the 
problem completely, it prevents one from going to small lattice spacings. In our formu- 
lation we do not have that constraint on the lattice spacing. However, the price we pay 
is that we resolve the problem of the linearization partially, not completely. Thus, it 
is of great importance for one to see to what extent we've resolved the problem and, if 
possible and necessary, to find a way to get around it completely. Up to this point we 
are not able to answer this question, but this is one of the crucial steps we should take 
towards a more sensible supersymmetric model on a lattice within our formulation. 

B J\f=l super-algebra in Majorana representation 

In this appendix we give the super-algebra in Majorana representation in D + 1 dimen- 
sional light-cone coordinates where D = 1,2,3. 

In Majorana representation Majorana spinors have real component fields, and can 
be written as 



where 6l, 9 r are left-moving, right-moving spinors with real components. This implies 
that the supercharge Q is also a Majorana spinor with real components of the form 
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where the integration is taken over the D spatial dimensions, is the supercurrent, 
which is a Majorana spinor. 

In terms of the Majorana super-charge, the super-algebra is given by 

{Q,Q} = 2T»P„ (22) 

where Q = Q^T° in any representation, and thus Q = Q T T° in Majorana representation. 

B.l D=l 

For 1+1 dimensional case, we have T° = a 2 and T 1 = ia 1 , so that 



10 -l Fl . ( 

% 



r o + r i 



and Q = i(Q~ , -Q+). Thus, Eq. (22) reads 



r - = 



r o_ r i 



^2 



or 



{Q ± , Q ± } = 2v / 2P ± , {q + , g-} = o. 



-V2 




-2y/2P* 




B.2 D=2 

In this case T° = <r 2 , T 1 = ia 1 and r 2 = Y 1 - = ia 3 . Therefore, 



or 



{Q ± ,Q ± } = 2V2P ± , {Q + ,Q'} = -2P A 



( -2P L -2 V / 2P H 



2P J 



B.3 D=3 

In 3+1 dimensions Majorana spinors have four components and thus the supercharge 
can be written as 



Q 



Q" 
Q- 



( Qt \ 
Qt 
Qi 
V Q2 J 



Q = Q T Y° = i(Q 2 , -Ql, Qt, -Qt) 



Gamma matrices are 4x4 matrices given by 



r+ 



r° + r 3 



o \ r°-r 3 

V2a 2 ) ' - v /2 



r 3 = 

^2 




-(72 
<7 2 
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Then Eq. (22) yields 



{Q,Q} = 2T»P fl = 2 



2i 



( -P 2 
-P 1 


\ V2P- 



ia 1 P 1 + ia 3 P 2 V2a 2 P_ 
v4j 2 P + i^Pi + i(r 3 P 2 
-V2P+ \ 



-P 1 

P 2 V2P- 

-V2P- -P 2 

-P 1 





-P 1 
P 2 



Hence, we find 



{Qt,Q^} = 2V2P ± 5 aP , 

{Qt, Qi) = -{Qt, Q2} = 2P 1 , {Qt, Q2} = {Qt, Qi) = -2P 2 , 

where a, (3 = 1,2. Note that if P 1 = P 2 = 0, then this algebra coincides with the one 
for M=2 SUSY in 1+1 dimensions also known as W=(2,2) SUSY. 
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